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Abstract 

An update of the search for sleptons, neutralinos and charginos in the context 
of scenarios where the lightest supersymmetric particle is the gravitino and the 
next-to-lightest supersymmetric particle is a slepton, is presented, together with 
the update of the search for heavy stable charged particles in light gravitino 
scenarios and Minimal Supersymmetric Standard Models. Data collected in 
1999 with the DELPHI detector at centre-of-mass energies around 192, 196, 
200 and 202 GeV were analysed. No evidence for the production of these 
supersymmetric particles was found. Hence, new mass limits were derived at 
95% confidence level. 
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1 Introduction 

In 1999, the centre-of-mass energies reached by LEP ranged from 192 GeV to 202 
GeV, and the DELPHI experiment collected an integrated luminosity of 228.2 pb _1 . 
These data were analysed to update the searches for sleptons, neutralinos and charginos 
in the context of Gauge Mediated Supersymmetry Breaking (GMSB) models [[TJQ already 
performed at lower energies. 

In these models the gravitino, G, is the lightest supersymmetric particle (LSP) and 
the next-to-lightest supersymmetric particle (NLSP) can be either the neutralino, x?> or 
the sleptons, I |||3],[|J|[7|| . The data were analysed under the assumption that the NLSP 
is a slepton. Depending on the magnitude of the mixing between the left and right gauge 
eigenstates, tr and tl, there are two possible scenarios. If the mixing is largeQ, fx (the 
lighter mass eigenstate) is the NLSP. However, if the mixing is negligible, f\ is mainly 
right-handed || and almost mass degenerate with the other sleptons. In this case, the 
e.R and (Ir three body decay (I— > fir/ with f\—> r G), is very suppressed, and e# and 
jiR decay directly into IG. This scenario is called sleptons co-NLSP. 

Due to the coupling of the NLSP to G, the mean decay length, L, of the NLSP can 
range from micrometres to metres depending on the mass of the gravitino |J (mg): 

L = 1.76 x 10-\l(^-) 2 - 1 ( m > VV-^* Y cm (1) 
V W UOOGeV/cV VleV/cV V; 

For example, for ttiq < 250 eV/c 2 , or equivalently, for a SUSY breaking scale of 
\J~F~ < 1000 TeV (since both parameters are related ||10|), the decay of a NLSP with mass 
greater than for example 60 GeV/c 2 can take place within the detector. This range of 



'F is in fact consistent with astrophysical and cosmological considerations |11| , |12 • 
In this work the results of the searches reported in are updated and the search 
for charginos is extended to the sleptons co-NLSP scenario. Moreover, the update of the 



search for heavy stable charged particles presented in JTj] is also performed. Heavy stable 
charged particles are predicted not only in GMSB models but also in Minimal Supersym- 
metric Standard Models (MSSM) with a very small amount of R-parity violation, or with 
R-parity conservation if the mass difference between the LSP and the NLSP becomes 
very small. In these models the LSP can be a charged slepton or a squark and decay 



with a long lifetime into Standard Model particles [[TIJ. Therefore, updated lower mass 



limits on heavy stable charged particles, under the assumption that the LSP is a charged 
slepton, will be provided in this letter within both models, GMSB and MSSM. 
The first search looks for the production of Xi pahs in the f% NLSP scenario 

e + e" - XiXi - r+T-f+T- - t + Gt-t + Gt- (2) 

and in the sleptons co-NLSP scenario 

e + e~ -> XiXi -> Ir~1~IrI'~ -> l + Gri' + Gl'- (3) 

with I = e, /i, r and BR(x ( j ) — > 11) = 1/3 for each leptonic flavour. In the former case, 
neutralino pair production would mainly lead to a final state with four tau leptons and 
two gravitinos, while in the case of a co-NLSP scenario, the final signature would contain 
two pairs of leptons with possibly different flavour and two gravitinos. 

tin GMSB models large mixing occurs generally in regions of tan /3 > 10 (the ratio of the vacuum expectation values of 
the two Higgs doublets) or \fi\ > 500 GeV/c 2 (/x is the Higgs mass parameter). 



2 



The second search concerns / pair production followed by the decay of each slepton 
into a lepton and a gravitino: 

e + e- -> l + r -> l + Gl~~G (4) 

This search has been performed within the fi NLSP scenario (I = fi) and the sleptons 
co-NLSP scenario (I = In). The signature of these events will depend on the mean decay 
length of the NLSP, or equivalently, on the gravitino mass. Therefore, if the decay length 
is too short (1 eV/c 2 < < 10 eV/c 2 ) to allow the reconstruction of the slepton, 
only the corresponding lepton or its decay products will be seen in the detector, and the 
search will then be based on the track impact parameter. If the slepton decays inside the 
tracking devices (10 eV/c 2 < < 1000 eV/c 2 ), the signature will be at least one track 
of a charged particle with a kink or a decay vertex. However, for very heavy gravitinos 
( m G ~ 1000 eV/c 2 ), the decay length is large and the slepton decays outside the detector. 
The pair production of such long-lived or stable particles yields a characteristic signature 
with typically two back-to-back charged heavy objects in the detector. Finally, for very 
light gravitino masses (m^ < 1 eV/c 2 ), the decay takes place in the primary vertex and 
the results from the search for sleptons in gravity mediated (MSUGRA) models can be 
applied fl6|| . 



In the parameter space where the sleptons are the NLSP, there are specific regions 
where the chargino is light enough to be produced at LEP ||. Therefore, the third 
search looks for the pair production of lightest charginos in the f\ NLSP scenario 

e + e~ — > XiXi ~^ T^v T f^V T — > t + Gv t t~GV t (5) 

and sleptons co-NLSP scenario 

e + e" -> xtXi -> l%yil r R V V -> 1+GvJ-GVv. (6) 

The analysis is divided into four topologies according to the mean lifetime of the slepton as 
explained in the previous paragraph: two acoplanar leptons with respect to the beam pipe 
with missing energy (MSUGRA models), at least one track with large impact parameter 
or a kink, or at least one track corresponding to a very massive stable charged particle. 

The data samples are described in section ^. The efficiencies of the different selection 
criteria and the number of events selected in data and in the Standard Model background 
are reported in Finally, the results are presented in section [|. 



2 Data sample and event generators 

All searches are based on data collected with the DELPHI detector during 1999 at 
centre-of-mass energies around 192, 196, 200 and 202 GeV. The total integrated lumino- 



sity was 228.2 pb . A detailed description of the DELPHI detector can be found in | i7| 
and the detector performance in Jl 



To evaluate the signal efficiencies and background contamination, events were genera- 



ted using different programs, all relying on JETSET 7.4 ||19|| , tuned to LEP 1 data for 
quark fragmentation. 

The program SUSYGEN [^l|] was used to generate neutralino pair events and their sub- 
sequent decay products. In order to compute detection efficiencies, a total of 90000 events 
were generated with masses 67 GeV/c 2 < m ?1 + 2 GeV/c 2 < m^o < \fsj1 and at the four 
centre-of-mass energies. 
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Slepton pair samples of 99000 and 76500 events at 196 GeV and 202 GeV centre-of- 
mass energies respectively were produced with PYTHIA 5.7| [[19] with staus having a mean 
decay length from 0.25 to 200 cm and masses from 40 to 100 GeV/c 2 . Other samples of f 
pairs were produced with SUSYGEN for the small impact parameter search with rrif from 
40 GeV/c 2 to 100 GeV/c 2 . 

For the search for heavy stable charged particles, signal efficiencies were estimated on 
the basis of about 50000 simulated events. Pair produced heavy smuons were generated 
at energies between 192 GeV and 202 GeV with SUSYGEN, and passed through the de- 
tector simulation as heavy muons. The efficiencies were estimated for masses between 
10 GeV/c 2 and 97.5 GeV/c 2 . 1000 events were generated per mass point. 

SUSYGEN was also used to generate Xi P a i r production samples and their decays 
at 192 GeV and 202 GeV centre-of-mass energies. In order to compute detection effi- 
ciencies, a total of 45 samples with 500 events each were generated with rriQ at 1, 100 
and 1000 eV/c 2 , m fl + 0.3 GeV/c 2 < m~+ < y/s/2 and m fl > 65 GeV/c 2 . Samples with 



smaller Am = m-+ — were not generated because in this region the xt does not 



'x 



decay mainly to f\ and v T but into W and G. 

The background process e + e~ — > qq(n , ~f) was generated with PYTHIA 6.125, while 
K0RALZ 4.2 [H] was used for fi + fi~(j) and r + r"(7). The generator BHWIDE [||] was 



used for e + e — > e + e events. 

Processes leading to four-fermion final states were generated using EXCALIBUR 



1.08 11 and GRC4F 25 



Two-photon interactions leading to hadronic final states were generated using 
TW0GAM |f26[ , including the VDM, QPM and QCD components. The generators of Berends, 



Daverveldt and Kleiss |27| were used for the leptonic final states. 

The cosmic radiation background was studied using the data collected before the 
beginning of the 1998 LEP run. 

The generated signal and background events were passed through the detailed simu- 
lation |18] of the DELPHI detector and then processed with the same reconstruction and 
analysis programs used for real data. 



3 Data selection 

3.1 Neutralino pair production 

The selection criteria used in the search for neutralino pair production in the fx NLSP 



scenario and in the sleptons co-NLSP scenario, were described in detail in fl3l , |28|j . The 
main two differences between these two cases are that the mean number of neutrinos 
carrying away undetected energy and momentum and the number of charged particles 
per event are considerably bigger for the fi NLSP scenario. 

After applying the selection criteria to the search for these topologies, six events passed 
the search for neutralino pair production in the f\ NLSP scenario, with 3.36±0.98 Stan- 
dard Model (SM) background events expected. Four events passed the search for neu- 
tralino pair production in the sleptons co-NLSP scenario, with 4.39±0.51 SM background 
events expected. Efficiencies between 20% and 44% were obtained for the signal events. 



t Another two samples of 1000 events with rrif = 60 GeV/c 2 and mean decay lengths of 5 and 50 cm, were generated 
using SUSYGEN to cross check with the PYTHIA results. The same efficiencies were found within a ±2% difference. 
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3.2 Slepton pair production 

This section describes the update of the search for slepton pair production as a function 
of the mean decay length. The details of the selection criteria used to search for the 
topologies obtained when the NLSP decays inside the detector volume were described 
in |T3|j28|j29|| . Likewise, the selection criteria used to search for heavy stable charged 
particles were described in detail in |T^ , |30| . The efficiencies were derived for different 



I masses and decay lengths by applying the same selection to the simulated signal events. 
3.2.1 Search for secondary vertices 



This analysis exploits a feature of the I — > IG topology when the slepton decays inside 
the tracking devices, namely, one or two tracks originating from the interaction point 
and at least one of them with a secondary vertex or a kink. After applying the selection 
criteria to search for this topology, two events in real data were found to satisfy all 
the requirements, while 0. 79^012 were expected from SM backgrounds. One event was 
compatible with a 77 — > r + r~ with a hadronic interaction in the Inner detector. The 
other one was compatible with a e + e~ — > t + t~ event where one of the electrons (a 
decay product of the r), after radiating a photon, was reconstructed as two independent 
tracks. 

The secondary vertex reconstruction procedure was sensitive to radial decay lengths, 
R, between 20 cm and 90 cm. The Vertex detector and the Inner detector were needed to 
reconstruct the f and the Time Projection Chamber to reconstruct the decay products. 
The shape of the efficiency distribution was essentially flat as a function of R decreasing 
when the f decayed near the outer surface of the Time Projection Chamber. The de- 
crease was due to inefficiencies in the reconstruction of the tracks coming from the decay 
products of the r. The search for events with secondary vertices had an efficiency of 
(52.8±2.0)% for f with masses between 40 and 100 GeV/c 2 , with a mean decay length 
of 50 cm. 

The same selection criteria were applied to smuons and selectrons. The efficiency was 
(56.5±2.0)% for m~ m between 40 to 100 GeV/c 2 , and for selectrons it was (38.1±2.0)% 
in the same range of masses. The efficiency for selectrons was lower than for staus or 
smuons due to an upper cut on total electromagnetic energy at the preselection level. 

3.2.2 Large impact parameter search 



To investigate the region of lower gravitino masses the previous search was extended to 
the case of sleptons with mean decay length between 0.25 cm and approximately 10 cm. 
In this case the I is not reconstructed and only the / (or the decay products in the case of 
f ) is detected. The impact parameter search was only applied to those events accepted 
by the same general requirements as in the search for secondary vertices, and not selected 
by the secondary vertex analysis. 

The maximum efficiency was (31.0±2.0)% corresponding to a mean decay length of 
2.5 cm. The efficiency decreased sharply for lower decay lengths due to the requirement on 
minimum impact parameter. For longer decay lengths, the appearance of reconstructed 
f in combination with the cut on the maximum number of charged particles in the 
event caused the efficiency to decrease smoothly. This decrease is compensated by a 
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rising efficiency in the search for secondary vertices. For masses above 40 GeV/c 2 no 
dependence on the f mass was found far from the kinematic limit. 

The same selection was applied to smuons and selectrons. For smuons the efficiency in- 
creased to (60.0±2.0) % for a mean decay length of 2.5 cm and masses over 40 GeV/ c 2 since 
the smuon always has a one-prong decay. For selectrons the efficiency was (36.3±2.0)% 
for the same mean decay length and range of masses. 

Trigger efficiencies were studied simulating the DELPHI trigger response to the events 
selected by the vertex search and by the large impact parameter analysis, and were found 
to be around 99%. 

Two events in the real data sample were selected, while 1.77^o2i were expected from 
SM backgrounds. Both events are compatible with e + e~ — > t + t~ events where one of 
the electrons (a decay product of the r), after radiating a photon, was reconstructed only 
by the Inner detector and Time Projection Chamber detectors giving a very large impact 
parameter track. 

3.2.3 Small impact parameter search 



The large impact parameter search can be extended further to mean decay lengths 
below 0.1 cm. Here only the main points of the analysis and some changes with respect 
to previous ones are recalled. In low multiplicity events two hemispheres were defined 
using the thrust axis. The highest momentum, good quality (Ap/p < 50%), particle 
tracks in each hemisphere were labelled leading tracks. The impact parameters from the 
beam spot, bi and b 2 , of the leading tracks in the R(p plane were used to discriminate 
against SM backgrounds. The same selection criteria described in references [0,^1 were 



applied. However, some extra selection was added in order to reduce the background 
from detector noise or failure. 

In order to preserve the efficiency in the region of decay length > 10 cm, where the 
f can be observed as a particle coming from the primary vertex and badly measured 
owing to its limited length, further requirements on the track quality were applied only 
to the leading track with the larger impact parameter. This particle was required to have 
a relative momentum error < 30% and the track to be measured at least either in the 
Time Projection Chamber or in all of the other three track detectors in the barrel (Vertex 
detector, Inner detector and Outer detector). 

The efficiency of the search did not show any significant dependence on the f mass for 
masses over 40 GeV/c 2 and it could be parameterized as a function of the f decay length 
in the laboratory system. The maximum efficiency was ~ 38% for a mean decay length 
of ~ 2 cm; the efficiency dropped at small decay lengths (~ 10% at 0.6 mm). 

The same selection criteria were used to search for smuons as reported in |T3 



The maximum efficiency reached for the smuon search was ~ 43% at ~ 2 cm 
mean decay length. To search for selectrons, in order to increase efficiency, the cut 
(Ei + E 2 ) < 0.7 Ef, eam (where E±, E 2 are the electromagnetic energy deposits associated 
to the leading tracks) was not applied. The Bhabha events that survived the selection, 
when the previous rejection cut was not applied, were those where at least one of the 
electrons underwent a secondary interaction, thus acquiring a large impact parameter. 
However, it was found that in these cases the measured momentum of the electron was 
smaller than the electromagnetic energy deposition around the electron track. Therefore, 
the cut (Ei/pi + E 2 /p 2 ) < 2.2 was used for the selectron search. The maximum efficiency 
reached in the selectron search was ~ 35% at ~ 2 cm mean decay length. 



6 



Requiring y/b\ + b\ > 600 /an, the number of events selected in the data was 5 in the 
f and fi search, and 4 in the e search, while 5.05±0.39 events were expected from the SM 
background in both searches. All of the selected candidates were compatible with SM 
events. 



3.2.4 Heavy stable charged particle search 

The analysis described in |L4j] has been applied for each of the four centre-of-mass 
energies 192, 196, 200 and 202 GeV. A careful run selection ensured that the Ring Imaging 
CHerenkov (RICH) detectors were fully operational because the method used to identify 
heavy stable particles relies on the lack of Cherenkov radiation in DELPHI'S RICH de- 
tectors. The luminosities analysed after the run selection are shown in Table |l|. Signal 
efficiencies were estimated from Monte Carlo by simulating heavy sleptons with SUSYGEN 
and passing them through the detector simulation as heavy muons. The background 
given in Table [l] was estimated from data itself by counting the number of tracks passing 
the individual selection criteria. Only events with two or three charged particles were 
considered. Events were selected, if they contained at least one charged particle with: 

1. a momentum above 5 GeV/c, a high ionization loss and no photons in the gas 
radiator were associated to the particle (gas veto) or, 

2. a momentum above 15 GeV/c, an ionization loss at least 0.3 below the expectation 
for a proton and surviving the gas veto or, 

3. a momentum above 15 GeV/c, surviving the gas and the liquid RICH veto. 

An event was also selected if both event hemispheres contained particles with both either 
a high ionization loss or a gas veto, or both having a low ionization loss. 

No candidate events were selected in data. As an example, Figure |T] show the data and 
the three main search windows for the search at an energy of 202 GeV. The expectation 
for a 90 GeV/c 2 mass signal is also shown. For particle masses below 60 GeV/c 2 the signal 
efficiencies are of the order of 30%, and rise with increasing mass to about 78%. Then the 
efficiency drops when approaching the kinematical limit due to saturation effects, and it 
is assumed to be zero at the kinematical limit. 





192 GeV 


196 GeV 


200 GeV 


202 GeV 


C (pb- 1 ) 


26.3 


69.7 


87.1 


40.4 


background 


0.12±0.04 


0.18±0.04 


0.31±0.06 


0.1±0.03 


observed 















Table 1: Luminosities analysed and selected events for each of the four centre-of-mass 
energies in the search for heavy stable charged particles. 



3.3 Chargino pair production 

The search for the lightest chargino depends on the slepton lifetime, or equivalently 
on the gravitino mass as already stated. For vtiq < 1 eV/c 2 , Xi decays at the vertex 
and the final state is two acoplanar leptons with missing energy. In this case the search 
for charginos and for sleptons in gravity mediated scenarios (MSUGRA) can be applied. 
The details of the search for charginos in MSUGRA models can be found in [[H]] . The 
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efficiencies obtained using these analyses were 13-36% for the fi NLSP scenario and 15- 
29% for the sleptons co-NLSP scenario. The number of selected events in data was 39, 
while the background expected from SM was 37.581q 90 f° r the f% NLSP scenario. For the 
sleptons co-NLSP scenario the number of candidates in data was 81, while the number 
of events expected from SM was 79.7±3.9. For uiq between 1 eV/c 2 and 1000 eV/c 2 , 
xf has intermediate mean decay lengths and the final topologies are events with kinks or 
large impact parameter tracks. The efficiencies obtained with these analyses were 25-56% 
in the f\ NLSP scenario, and 41-56% in the sleptons co-NLSP scenario. Four events in 
real data were found to satisfy all the conditions required in the search for charginos 
with intermediate mean decay lengths, while 2.56^Q24 events were expected from SM 
backgrounds. These results apply to both scenarios f\ NLSP and sleptons co-NLSP, 
since the same selection criteria were applied to search for staus, smuons and selectrons. 

Finally, for ttiq > 1000 eV/c 2 the event topology is at least one track corresponding to 
a very massive stable charged particle. Therefore the search for stable charged particles 
was applied. In this case the efficiency is mainly affected by the momentum of the slepton 
because the method used to identify heavy stable particles relies on the lack of Cherenkov 
radiation in DELPHI'S RICH detectors as already stated in section |3 . 2 . 4| . To remove SM 
backgrounds, low momentum particles are removed, thus reducing the efficiency for higher 
chargino masses, especially in the region where the mass difference (Am) between the 
NLSP and the LSP is small. Therefore for this analysis the efficiencies vary from 0% for 
Am = 300 MeV/c 2 to 62% for Am = 20 GeV/c 2 . No candidates in the data passed the 
selection cuts, while 0.71±0.09 events were expected from background. 



4 Results and interpretation 

Since there was no evidence for a signal above the expected background, the number 
of candidates in data and the expected number of background events were used to set 
limits at the 95% confidence level (CL) on the pair production cross-section and masses 
of the sparticles searched for. The model described in reference [|J was used to derive 
limits within the GMSB scenarios. This model assumes radiatively broken electroweak 
symmetry and null trilinear couplings at the messenger scale. The corresponding para- 
meter space was scanned as follows: 1 < n < 4, 5 TeV < A < 90 TeV, 1.1 < M/A < 10 9 , 
1.1 < tan/5 < 50, and sign(fi) = ±1, where n is the number of messenger generations 
in the model, A is the ratio between the vacuum expectation values of the auxiliary 
component and the scalar component of the superfield and M is the messenger mass 
scale. The parameters tan j3 and [i are defined as for MSUGRA. The limits presented 
here are at y/s = 202 GeV after combining the results of the searches at lower centre-of- 



mass energies with the likelihood ratio method [32]. 



4.1 Neutralino pair production 

Limits for neutralino pair production cross-section were derived in the f\ NLSP 
and sleptons co-NLSP scenarios for each (m^o,m z ~) combination. For the f\ NLSP 
case, the combination took into account the results from the LEP runs from 1996 (for 
\fs > 161 GeV) to 1999. The limits for the production cross-section allowed some sec- 
tors of the (m^o , mj) space to be excluded. In order to exclude as much as possible of the 
mass plane, the results from two other analyses were taken into account. The first is the 
search for slepton pair production in the context of MSUGRA models. In the case where 
the MSUGRA Xi is massless, the kinematics correspond to the case of I decaying into a 
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lepton and a gravitino. The second is the search for lightest neutralino pair production 
in the region of the mass space where Xi is the NLSP j33| (the region above the diagonal 
line in Figure ^], i.e. > m^o). Within this zone, the neutralino decays into a gravitino 
and a photon. 

As an illustration, Figure |2| presents the 95% CL excluded areas for rriQ < 1 eV/c 2 in 
the m^o vs m fl plane for the fx NLSP. The negative-slope dashed area is excluded by the 

analysis searching for neutralino pair production followed by the decay Xi ~^ The 
point-hatched area is excluded by the direct search for slepton pair production within 
MSUGRA scenarios. 



4.2 Slepton pair production 

The results of the search for slepton pair production are presented in the (mQ,mj) 
plane combining the two impact parameter searches, the secondary vertex analysis and 
the stable heavy lepton search, and using all DELPHI data from 130 GeV to 202 GeV 
centre-of-mass energies. 

The fi pair production cross-section depends on the mixing in the stau sector. There- 
fore, in order to put limits on the fx mass the mixing angle had to be fixed. The results 
presented here correspond to the case when there is no mixing between the tr and fx, 
thus fi is a pure right-handed state (Figure |3|-a). In the case which corresponds to 
a mixing angle which gives the minimum fi pair production cross-section and at the 
same time maintains m~ > 0, the given limit was reduced by ~ 1 GeV. Therefore, 
within the f\ NLSP scenario, the impact parameter and secondary vertex analyses ex- 
tended the limit m fR > 75 GeV/c 2 for m d < 1 eV/c 2 , set by MSUGRA searches pj, 
up to ttiq = 400 eV/c 2 , reaching the maximum excluded value of vnf R = 92 GeV/c 2 for 
mg= 200 eV/c 2 . For tjIq > 250 eV/c 2 the best lower mass limit was set by the stable 
heavy lepton search. 

Within the sleptons co-NLSP scenario, the cross-section limits were used to derive 
lower limits for (Figure ^|-b) and e# (Figure ||-c) masses at 95% CL. The JIr pair pro- 
duction cross-section is model independent, however, the e# pair production cross-section 
is a function of the GMSB parameters due to the exchange of a Xi m the t-channel. 
Therefore, in order to put limits on the cr mass, the aforementioned region of the GMSB 
parameter space was scanned and, for each selectron mass, the smallest theoretical pro- 
duction cross-section was chosen for comparison with the experimental limits. For gra- 
vitino masses below a few eV/c 2 , the experimental limits are the ones corresponding to 
the search for selectrons in MSUGRA models. 

Assuming mass degeneracy between the sleptons, (Figure 0-d), these searches extended 
the limit l R > 80 GeV/c 2 set by MSUGRA searches [0] for very short NLSP lifetimes, 
up to vn,Q= 700 eV/c 2 . For the MSUGRA case no lepton combination exists, so the best 
limit from the JIr has been used. The maximum excluded value of ^i R = 94 GeV/c 2 was 
achieved for vn,Q= 200 eV/c 2 . For uiq > 700 eV/c 2 the best lower mass limit was set 
by the stable heavy lepton search. Ir masses below 35 GeV/c 2 were excluded by LEP 1 



data |j3~4| . In the case of Ir degeneracy, this limit improved to 41 GeV/c 



4.3 Chargino pair production 



The limits on the chargino pair production cross-section were used to exclude areas 
within the (m^ , tuj) plane for different domains of the gravitino mass combining results 
from all the centre-of-mass energies from 183 GeV to 202 GeV. Figure |] shows the 
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regions excluded at 95% CL in the (m~+,m fl ) plane (a) and (m-+,m^) plane (b). The 
positive-slope area is excluded for all gravitino masses. The negative-slope area is only 
excluded for ttiq > 100 eV/c 2 and the brick area for mg > 1 keV/c 2 . The areas below 
mf,= 73 GeV/c 2 (Figure |] (a)) and below mr = 80 GeV/c 2 (Figure |] (b)), are excluded 
by the direct search for slepton pair production in MSUGRA models ||16|| . The area of 
Am < 0.3 GeV/c 2 is not excluded because in this region the charginos do not decay 
mainly to fi and u T , but to W and G. Thus, if Am > 0.3 GeV/c 2 , the chargino mass 
limits are 95.2, 96.8 and 99 GeV/c 2 for = 1, 100 and 1000 eV/c 2 respectively, in 
the fi NLSP scenario. In the sleptons co-NLSP scenario the limits are 95.2, 98.6 and 
98.6 GeV/c 2 for = 1, 100 and 1000 eV/c 2 respectively. The limit at m^ = 1 eV/c 2 
is also valid for smaller masses of the gravitino, because they lead to the same final 
state topologies. The same argument is true for m^ > 1 keV/c 2 . The chargino mass 
limit decreases with decreasing because in scenarios with gravitino LSP, small stau 
masses correspond to small sneutrino masses (both are proportional to A) and hence to 
smaller production cross-sections due to the destructive interference between the s- and 
t-channels. It should be noticed that within the parameter space covered by this work, the 
lightest chargino is at least 40% heavier than the lightest neutralino. Thus, for gravitino 
masses up to ~ 1 eV/c 2 the search for neutralinos implies a model dependent lower limit 
on the lightest chargino of 125 GeV/c 2 . However, neutralinos were not directly searched 
for in heavier gravitino mass regions, therefore, a model dependent lower limit cannot be 
set in this case. Thus, the experimental lower limit of 98.6 GeV/c 2 remains valid. 

4.4 Heavy stable charged particle pair production 

The results presented in section [3.2.4| were combined with previous DELPHI results 
in this channel, and cross-section limits were derived as indicated in Figure [5[ From the 
intersection points with the predicted cross-sections for smuon or staus in the MSSM, 
left(right) handed smuons and staus can be excluded up to masses of 94.0(93.5) GeV/c 2 
at 95%CL. No limits are given on selectrons here, because the cross-section can be highly 
suppressed by an additional t-channel sneutrino exchange contribution. 

4.5 Limits on the GMSB parameter space 

Finally, all these results can be combined to produce exclusion plots within the 
(tan p, A) space. The corresponding parameter space was scanned as follows: 1 < n < 4, 
5 TeV < A < 90 TeV, 1.1 < M/A < 10 9 , 1.1 < tan/? < 50, and sign(n) = ±1. As an 
example, Figure |6] shows the zones excluded for n =1 to 4 for m^ < 1 eV/c 2 , which cor- 
responds to the NLSP decaying at the main vertex. The shaded areas are excluded. The 
areas below the dashed lines contain points of the GMSB parameter space with Xi NLSP. 
The areas to the right (above for n — 1) of the dashed-dotted lines contain points of the 
GMSB parameter space where sleptons are the NLSP. It can be seen that the region of 
slepton NLSP increases with n. The contrary occurs to the region of neutralino NLSP. 
A limit could be set for the variable A at 17.5 TeV. 

5 Summary 

Lightest neutralino, slepton and chargino pair production were searched for in the 
context of light gravitino models. Two possibilities were explored: the f\ NLSP and the 
sleptons co-NLSP scenarios. No evidence for signal production was found. Hence, the 
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DELPHI collaboration sets lower limits at 95% CL for the mass of the \i a ^ 86 GeV/c 2 
if rriQ < 1 eV/c 2 , and lower mass limits for the sleptons in all the gravitino mass range. 
The limit on the chargino mass is 95.2 GeV/c 2 for all in both scenarios, fi NLSP and 
the sleptons co-NLSP. 

Finally, mass limits for heavy stable charged particles were also derived within the 
MSSM. For these particles the DELPHI collaboration sets lower mass limits at 95% CL 
for the left (right) handed sleptons at 94.0 (93.5) GeV/c 2 . 
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Figure 1: (a) Normalised energy loss as a function of the momentum after the gas 
veto for the 202 GeV data, (b) Measured Cherenkov angle in the liquid radiator as a 
function of the momentum after the gas veto: if four photons or less were observed in 
the liquid radiator, the Cherenkov angle was set equal to zero. The rectangular areas 
in (a) indicate selections (1) and (2), and that in (b) shows selection (3). The selection 
criteria are explained in the text. Open circles are data. The small filled circles indicate 
the expectation for a 90 GeV/c 2 mass signal with charge ±e, resulting in a large dE/dx 
(upper plot) and no photons (except for a few accidental rings) in the liquid Cherenkov 
counter (lower plot). The solid lines with a mass signal value indicate the expectation 
for heavy stable sleptons. 
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Figure 2: Areas excluded at 95% CL with rriQ < 1 eV/c 2 in the m^o vs plane 
for n =1 to 4, using all data from 161 GeV to 202 GeV centre-of-mass energies. The 
positive-slope dashed area is excluded by this analysis. The negative-slope dashed area 
is excluded by the search for Xi ~ * jG, and the point-hatched area by the direct search 
for stau pair production in the MSUGRA framework. 
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Figure 3: Exclusion regions in the (mQ,m fR ) (a), (m^ra^) (b), (mQ,mg R ) (c) (mg 
(d) planes at 95% CL for the present analyses combined with the Stable Heavy Lepton 
(SHL) search and the search for I in gravity mediated models (MSUGRA), using all 
DELPHI data from 130 GeV to 202 GeV centre-of-mass energies. The positive-slope 
hatched area shows the region excluded by the combination of the impact parameter and 
secondary vertex searches. The dashed line shows the expected limits. 
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Figure 4: Areas excluded at 95% CL in the (m-^m^J plane (a) and (m-+,m^) plane 
(b) using all DELPHI data from 183 GeV to 202 GeV centre-of-mass energies. The 
positive-slope area is excluded for all uiq. The negative-slope area is excluded only for 
m G > 100 eV/c 2 and the brick area for ttlq > 1 keV/c 2 . The grey area is excluded by 
the direct search for slepton pair production within MSUGRA models. 
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Figure 5: Predicted production cross-section for left and right handed stable smuons 
(staus) as a function of the particle mass. The cross-section limit indicated in the figure 
has been derived using all DELPHI data between 130 and 202 GeV. 
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Figure 6: Shaded areas in the (tan/3, A) plane are excluded at 95% CL. The areas below 
the dashed lines contain points of the GMSB parameter space with Xi NLSP. The areas 
to the right (above for n = 1) of the dashed-dotted lines contain points of the GMSB 
parameter space were sleptons are the NLSP. 



